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Abstract 

We find that studying the simplest of the coupled non-equilibrium 
growth equations of Barabasi by self-consistent mode coupling requires 
the use of dressed vertices. Using the vertex renormalization, we find a 
roughening exponent which already in the leading order is quite close 
to the numerical value. 



PACS number(s):05.10.Gg,05.40.-a,05.70.-a,64.60.Ht 



Models of interfacial growth have attracted a tremendous amount of at- 
tention since the pioneering work of Kardar, Parisi, Zhang (KPZ) [1,2]. A 
variety of interesting issues are associated with the KPZ equation and they 
have given rise to a variety of novel techniques [3]. Among the first analytic 
techniques used to tackle the KPZ system were the dynamic renormaliza- 
tion group (DRG) [4] and the self-consistent mode coupling scheme (SCMC) 
[5,6]. An important variant of the KPZ system was introduced by Ertaas 
and Kardar [7] and Barabasi [8] . This variant consisted of two coupled fields 
(as opposed to one field in KPZ) and is useful for studying the effects of 
a second non-equilibrium field on the growing interface. In these coupled 
field problems DRG has been employed, as also numerical techniques. One 
does not always get a stable fixed point with the DRG analysis which may 
sometimes indicate a failure of the perturbation scheme or may indicate a 
basic instability of the system. It is interesting to note that in many cases 
the exponents coming from the one-loop DRG analysis are not in very good 
agreement with the numerical analysis. This is exemplified in the simplest 
situation treated by Barabasi - an essentially linear system coupled according 
to 

dip j_, d 2 ip ^ 86 dip 
dt 2 dx 2 dx dx 2 
with < N lt2 (x 1 ,t 1 )N lt 2(x2,t 2 )>= 2D li2 5(x 1 -x 2 )6(t 1 -t 2 ). 



The field 6 satisfies the Edwards- Wilkinson equation and the field ip is 
coupled linearly via a gradient coupling to the 0-field. While the Edwards- 
Wilkinson model can be exactly solved, this is not true for eqn.(2) because of 
the multiplicative noise (note that 6 is a random field). The DRG recursion 
relations in this case yield for the roughening exponent a of the ^-field, the 
value a = 5/6 while the numerical value of a is nearly 0.68. The dynamical 
exponent z of the ^-field is found to be 2. Thus, in this case the dynamic 
exponent for both 6 and ip fields is found to be 2. We will call this "extended" 
dynamic scaling i.e. the time scale is independent of the nature of the field 
[9]. As it turns out, this is the only situation for this case. However this 
need not always be so. In another model considered by Ertaas and Kardar 
and Barabasi, 



there are two possibilities: 

i) z0 = = 3/2, this is the extended dynamic scaling and is found to be the 
correct situation for A > with Ai > 0, 

ii) = 3/2, but z^ = 2, this situation is obtained for A < with A x > and 
can be described as "weak" scaling [9]. For problems involving two or more 
coupled fields, one needs to differentiate between "extended" and "weak" 
scaling. 

In the one-dimensional KPZ, the perturbative DRG is exact (due to the 
existence of a fluctuation-dissipation relation), but this is not true for the 
coupled system in one dimension. The self-consistent mode coupling (SCMC) 
which has been reasonably succesful for the KPZ, has never been attempted 
in the coupled system. In this note, we apply the SCMC to the coupled sys- 
tem to see if it is a quantitatively better scheme than the perturbative DRG. 
In the process, we find something quite unusual. In all known situations, 
SCMC has been succesful in cases where the vertex is not renormalized. 
This, in contrast, is a situation where the momentum dependence of the 
dressed vertex is absolutely essential. This is what makes the application of 
SCMC interesting in this problem and should act as a prototype for situa- 
tions where dressed vertices are unavoidable. Writing eqns. (1) and (2) in 
momentum space, we have 

<f>(k) = -r x k 2 <f)(k) + Ni(k) (5) 

tp(k) = -r 2 k 2 i/j(k) - A P( k - P) <t>{p)^{k - P) + N 2 (k) (6) 

p 

with < N 1:2 (h, 0J\)N lt2 (k 2 , oj 2 ) > = 2£>i 2 5(k 1 + k 2 ) 5(oj 1 + uj 2 )A 



The basic elements of the calculation are the Green's functions G^k^u) 
and G^(k, u), the correlation functions C^k, uj) and C^(k, uj) and the vertex 
function A(k, q,k — q). Obviously and are exactly known and are given 
by 

G^(k,uj) = -iuj + T x k 2 (7) 

while for the ip-fie\d 

G^(k,uj) = -iuj + T 2 A; 2 + E{k,u) (9) 

on 

UJ A + 1 2 K 



and 



A(k,p,k-p) = A + A(k,p) (11) 
The self-energy T,(k,u) is found at the dressed one level to be given by 



T,(k,u) = X 2 [ < ^- < ^-kp 2 (k-p)C lj> {p,u')G 1 p{k-p,uj - u') 

X 2 D 1 f dp fc(fc-p) 

2^ -iu + T lP 2 + Yjjc -p) + T 2 (k - p) 2 



;i2) 



where we have used eq.(8) and eq.(9) in the Lorentzian approximation, 
i.e. during the frequency convolution, H^(k,u) has been replaced by its zero 
frequency form. 

Our first observation is that within the extended dynamic scaling, we 
expect — We need to examine if this is self-consistent. Setting E(fc) = 
Tk 2 , we have 

Yk 2 = f ^ Kk-p) 

T 1 J 2ir r lP 2 + f 2 (k-p) 2 1 1 

where f 2 = T 2 + Y. 
The long wave length property (k — > 0) of the integral on the right hand 
side is best seen by changing to the symmetric variables p' = + p which 

gives the o{k 2 ) contribution of the integral to be k 2 / % fftwi) V 1 ' 
This integral is divergent and needs to be cut-off at o(k), which spoils the 
k 2 behaviour. The only way this can be prevented is by setting 3I\ = T 2 , 
which makes the o(k 2 ) contribution of S vanish, i.e. implies T = and this 
establishes 

3r\ = r 2 (14) 

which is in exact agreement with the earlier work of Barabasi. 
We now discuss the correlation function. The diagram with bare vertex 
is shown in Fig. la and leads to 

p 2 {k - p) 2 C^p, u')C4k -p,u- u') (15) 
We now assume the scaling form 

C^) = ^f^/k 2 ) (16) 

which is consistent with the equal time correlation function, / ^C^(k, to) 
being k~ x ~ 2a . In the absence of A, a = 1/2 and the extra roughness pro- 
duced by this added noise is expected to raise a beyond 1/2. Our expectation, 



then is that the second term will dominate in eqn.(15). The power count of 
the second term in eqn.(15) shows that C^(k,uj) ~ k~ 4 ~ 2a which cannot 
match the power count of the left hand side for any value of a and hence 
a self-consistent formulation requires the vertex to acquire a momentum de- 
pendence. Dressing the vertex leads to the diagram in Fig. lb. Dropping the 
first term on the right hand side of eqn.(15) and dressing the vertex in the 
second leads to 



C f (k,u;) = \G f (k,uj) | 2 A J ~^-p 2 (k-p) 2 A(k,p,k-p) 

C^{p,u)')C^{k - p,w - u') (17) 

Since, we are interested in the k — > property of C^(/c, cu), the vertex 
that we need is Lim k ^ A(k,p, k — p) and if in this limit the vertex has the 
form A p where A is a constant, then the self-consistency in power counting 
is restored. The consistency of the amplitude is assured if (we evaluate the 
integral in eqn.(17) in the leading approximation [12] of k — > 0) 



1 = 



XDi r dp (k — p) 2 A(k,p,k — p) 



9IV J 2tt I k-p\ 1+2a {k-p) 2 + \p 2 
AAq-Di f dp 1 



9IV J 2ir \p 2a 
XA D 1 1 
121Y 2a - 1 



(18) 



We note in passing that the above momentum dependence of the vertex 
does not alter the conditions of eqn.(14). The self consistent equation for the 
vertex is shown in Fig. 2. Clearly 



A(h,q,k-q) = X J ^^p(p~ q)(k-p) 2 Gip&^G^p- q,cu) 

C^(k - p, uj)A(p, p-q, q)A(k -p,p-q,k-q) (19) 

Once again, the dressed vertex A that we are interested in corresponds to 
k — > 0. This vertex scales as q on the left hand side. Power count of the right 
hand side shows that it is a linear function of momentum as well and thus the 
two sides are matched in exponents. To impose the amplitude inconsistency, 
we evaluate the integral on the right hand side in the dominant region which 
corresponds to small values of p. This leads to 

1 = ^^ (20) 
Comparing with eqn.(13), we find 



tt 4+^ a59 (21) 

This is to be compared with the numerical value of a — 0.68. For a 
more careful analysis, eqns.(17) and (19) have to be solved numerically. This 
is an extremely formidable task because the dependence of A on the three 
variables (two independent) has to be charted out. 

As a final point, one would like to show that in this particular case, 
the weak scaling situation does not arise. If z^ were to be different from 
2, then for T,(k,u) to be at all relevant, has to be smaller than 2. This 
means eqn.(12), would at zero frequency become (we now include the vertex 
correction) 

- X2L>1 f dp A(fc,p)fc(fc-p) 

Hk) ~ TT J 2^ s(fc-p) { } 

Simple power counting shows that with A oc p, z^ = 2, which contradicts 
our starting assumption that < 2 and hence there is no self-consistent 
solution of the weak scaling variety. 

We have checked to ensure that for the extended scaling case, the self- 
consistent scheme does give the roughening exponent. Whether, the scheme 
can be made to work for the weak scaling situation is under consideration. 
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2 Figure Captions 

Fig. la The self-consistent equation for the correlator with bare vertex. The 
double thick line is the dressed correlator and the double straight line the 
propagator. The cross stands for the noise. 

Fig. lb The self-consistent equation for the correlator with dressed vertex. 
The double thick line is the dressed correlator and the double straight line 
the dressed propagator. The cross stands for the noise. 
Fig. 2 The self-consistent equation for the vertex. 



